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An analytical treatment has been developed to study some of the axisymmetric 
vortex breakdown and reconnection fluid dynamic processes underlying body- 
vortex interactions that are frequently manifested in rotorcraft and propeller- 
driven fixed-wing aircraft wakes. In particular, the presence of negative vorticity in 
the inner core of a vortex filament (one example of which is examined in this paper) 
subsequent to “cutting” by a solid body has a profound influence on the vortex 
reconnection, leading to analog flow behavior similar to vortex breakdown 
phenomena described in the literature. Initial vorticity distributions (three specific 
examples which are examined) without an inner core of negative vorticity do not 
exhibit vortex breakdown and instead manifest diffusion-like properties while 
undergoing vortex reconnection. Though this work focuses on laminar vortical 
flow, this work is anticipated to provide valuable insight into rotary-wing 
aerodynamics as well as other types of vortical flow phenomena. 


Nomenclature 

r Radial coordinate, origin at filament axis, m 

r* Nondimens ional radial coordinate, r*=r/s 

r Q Spherical finite-volume source terms’ “effective” initial radius, m 
r c0 “Finite core” vortex filament core size radius (at time equal to zero), m 
Re Vortex Reynolds number, Re = y/v 

s Axial distance of filament breakpoint (time equal zero) from origin, m 

t Time, sec 

t* Nondimens ional time parameter, t*=vt/s 2 

V Velocity vector, cylindrical coordinates, V = [v r v d v z ],m/sec 
V p Finite core vortex axial flow “wave front” propagation velocity, m/sec 

z Axial (along vortex filament axis) coordinate, origin at intersection of filament segments’ plane of 

symmetry and filament axis, m 
z* Nondimens ional axial coordinate, z*=z/s 

y Vortex filament initial circulation strength, m 2 /sec 

T Vortex circulation, m 2 /sec 

v Kinematic viscosity, m 2 /sec 

0 Angular coordinate, radians 

co Vorticity vector, co = [co r co 6 co z ], 1/sec 
ip Axisymmetric flow stream function 


1 Aerospace Engineer, Aeromechanics Branch, Flight Vehicle Research and Technology Division, Mail 

Stop 243-12, AIAA Associate Fellow. 



I. Introduction 


V ORTEX filament “cutting” by solid bodies occurs in a number of real world examples for rotary-wing 
aircraft. These examples include blade vortex interactions where vortices shed from one rotor blade 
collide with another blade, trailed tip vortex interactions/collisions with a helicopter airframe/tailboom in 
hover and low-speed flight, or, alternatively, proprotor trailed tip vortex collisions with wing surfaces for 
tiltrotor aircraft in airplane-mode (Fig. 1). The time (wake age) evolution of rotor or propeller wake trailed 
vortex filaments can have profound implications for rotary- wing performance, acoustics, and aerial vehicle 
aerodynamics. With the recent renewed emphasis on rotary- wing aircraft with multiple rotors (such as 
high-speed coaxial aircraft and quad-tiltrotors among other vehicle configurations), the fundamental 
understanding of body-vortex interaction fluid dynamics is becoming increasingly important. Though 
tremendous progress has been made with respect to theoretical and computational fluid dynamic studies, 
the full implications of unsteady vortical flow phenomena are still not well understood in the rotorcraft 
community. Therefore, body-vortex interaction is still very much an active area for research for rotorcraft. 
From an experimental perspective, only a limited amount of quantitative information exists in the literature. 
References 1-3, for example, provide qualitative discussion regarding the vortex filament “cutting” as 
manifested in the form of blade vortex interactions. References 4-7 provide similar qualitative discussion 
of vortex and fuselage/airframe/wing interactions. Additionally, a plethora of anecdotal evidence exists in 
form of observations and images of rotor wakes from various aircraft where the rotor helical vortices have 
been documented, usually by natural condensation in the air. But these images, though oftentimes visually 
stunning, e.g. Ref. 8, are not particularly useful in tackling the body- vortex interaction, as well as the larger 
rotor/airframe interactional aerodynamics, problem. 
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Fig. 1. Various Representative Body-Vortex Interactions; potentially including “cutting” of vortex 
filaments (from Ref. 9) 


Very little fundamental, including theoretical, work to date has been performed examining the post- 
cutting/collision unsteady fluid dynamics of vortex/solid body interactions. Most work has either 
concentrated on the vortex/body interaction process itself (see Ref. 10) or, alternatively, has looked at 
vortex-on-vortex interactions. A related topic is the study of ring vortices during vortex/vortex and 
vortex/body interactions; see Ref. 1 1, for example. Though flow visualization and measurement techniques 
are sufficiently maturing so as to examine in detail the rotary-wing/airframe vortex filament collision 
process - see, in particular, Refs. 4, 5, and 7 - it is clear that current rotor wake theoretical analysis is 
inadequate to accurately model the vortex breakdown or reconnection processes subsequent to body/vortex 
interaction. This work summarizes a new analysis approach to study these vortical flow processes, post- 
filament-cutting, during body/vortex interaction. 

Analytical descriptions of columnar vortex filament reconnection and breakdown phenomena — 
subsequent to “cutting” or, rather, collision with a solid body — were derived for four distinct and different 
initial vorticity distributions: instantaneous, uniform, parabolic, and dual-core (inner core with negative 
vorticity). Though all four initial vorticity distributions yielded several flow features in common with each 
other, there were also a number of flow phenomena that were unique to only a subset of the distributions 
studied. In particular, the dual-core model exhibited interesting vortical flow characteristics not only with 
regards to vortex filament reconnection but also manifested flow behavior analogous to that observed and 


(numerically) predicted for axisymmetric vortex breakdown. Figure 2 illustrates the three initial vorticity 
distributions studied. 



Fig. 2. Initial Vorticity Distributions Studied in Ref. 12: (a) uniform, (b) slightly parabolic, and (c) 
dual-core models 


Figure 3 illustrates the point that there is going to be varying degrees of vortex filament “breakage” (in 
form of filament distortion and separation of “breakpoints”) depending on the severity of the body/vortex 
interaction. Perpendicular collisions of vortex filaments with thin solid bodies will tend to minimize the 
separation distance between the filament segment breakpoints as well as minimize the filament distortion. 
Therefore, collision with thin solid bodies will tend to result in the “reconnection” of the vortex filaments. 
Conversely, interactions/collisions with bluff bodies will tend to result in large “breakpoint” separation 
distances and considerable filament segment distortion — therefore, most likely resulting in vortex 
“breakdown.” 
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Fig. 3. Hypothetical Vortex Filament Distortion/Breakage after Collision with Solid Bodies 


Undoubtedly, at this point of the discussion, while discussing “cut” or “broken” vortex filament 
segments, concern might be raised that the proposed analysis is perhaps in violation of Kelvin’s Theorem 
regarding the constancy of vortex circulation over (and the necessity of) closed paths. Or otherwise stated: 
“An important consequence of the (Kelvin) theorem is that a vortex-line cannot begin or end at any point in 
the interior of the fluid. Any vortex-lines which exist must either form closed curves, or else traverse the 
fluid, beginning and ending on its boundaries.” (Ref. 13). However, as has also been pointed out by many 



others (e.g. Ref. 14), Kelvin’s Theorem only holds absolutely true for inviscid flow. Further, the proposed 
analytical study of “cut” or “broken” vortex filament segments can be considered to fall into the same 
general category of gedanken (“thought”) investigation as the pioneering work of Klein’s Kaffeeloffel 
experiment (Ref. 15), wherein a plate set in motion in a fluid is “dissolved” so as to form and, 
subsequently, thereby leave behind a vortex sheet. This work will, of course, be shown to be consistent 
with all applicable governing fluid dynamic equations, as well as describing physically relevant vortical 
flow phenomena. 

Figure 4 illustrates an idealized version of the vortex filament in the intermediate process of vortex 
reconnection/breakdown, post-collision. Distortion of the filament as to deflection of the filament axis 
laterally from its undisturbed orientation is neglected. Initial axial flow (along the filament axis) is 
assumed to be zero as a result of the collision process: the filament, as it wraps around the solid body, will 
have to adhere to the no-flow boundary constraint. 



Fig. 4. Vortex Geometry Description 


Figure 5 illustrates the evolution of the vortex core radius iso-surfaces with time, based on the 
predictions for dual-core (and parabolic) vortex filaments discussed later in the paper. The core size 
distribution is one means by which the vortex reconnection can be observed. The increase in tangential 
velocity momentum in the intermediate region between the vortex breakpoints with increasing 
nondimensional time can be clearly seen in Fig. 5. Additionally, there is a localized increase in core size 

(relative to the “far- field” core size as z* -> oo) at the breakpoints. The temporal variation in core size in 
the immediate vicinity of the vortex breakpoints seems to grow, or fill-in, much faster for the parabolic 
initial vorticity distribution than for the dual- core distribution. Curiously, two of the initial vorticity 
distributions - the instantaneous and uniform-core distributions - do not exhibit a spatial (versus temporal) 
variation in vortex core size; instead in both these cases the vortex core contour is a constant diameter 
cylindrical surface. 




Table 1 presents a high level assessment of some of the flow phenomenology of the four different initial 
vorticity distributions studied. The hallmark flow features for vortex reconnection appear to be: the 
existence of axial flow “wave fronts” propagating away from the vortex filament “breakpoints” along the 
filament segment axes; “funnel- like” steam function contours in the outer flow (outside of the vortex core); 
and, consistent with the stream function predictions, the observation of spiral streakline, and particle trace, 
flow patterns also in the outer flow of the vortex filaments. These flow features are exhibited in the 
predictions made for all four initial vorticity distributions studied. Other noteworthy flow features that can 
be manifested as a part of the vortex reconnection process include: observation of complex and/or 
secondary flow structures in the vortex enstrophy contours at, or about, the vortex filament breakpoints 
early during the vortex reconnection process and spatial (versus temporal) variation of the vortex core 
radius. And, yet, these two additional flow features were not universally manifested by the initial vorticity 
distributions studied. Uniquely, the appearance of vortex breakdown phenomena, and secondary features 
in the circulation distribution in the intermediate region between vortex breakpoints, seems to only occur 
when there is negative vorticity in the inner core of the initial vorticity distribution (or, as yet to be 
established by future study, at least a region of reduced vorticity strength in the center of the vortex 
filament may be required). Finally, it should be noted that the analytical solutions for the finite-core 
vortices also provide valuable insight into unsteady columnar (i.e. Vq is only a function of r and t) vortex 
flow behavior for the special case of s = 0 ; in this regards, the work summarized in this paper is a valuable 
addition to the modest body of theoretical investigations into unsteady viscous vortex models. 




Table 1 - Examples of Vortex Reconnection Phenomenology 



Instantaneous 

Uniform-Core 

Parabolic-Core 

Dual-Core 

Axial “Wave Fronts” Propagating from Vortex 
Breakpoints 

X ( Vp — » oo) 

X 

X 

X 

“Funnel-Like” Stream Function Contours 
(Outer Flow) 

X 

X 

X 

X 

Spiral (Streakline and Particle -Trace) Flow 
Patterns (Outer Flow) 

X 

X 

X 

X 

Complex and/or Secondary Flow Structures in 
Enstrophy Contours 


X 

X 

(Occurs Earlier) 

X 

Spatial (versus Temporal) Variation of Vortex 
Core Radius 



X 

X 

Vortex Breakdown Bubbles 




X 

Secondary Features in Circulation “Diffusion” 
into Intermediate Region 




X 


The various types of flow phenomenon identified in Table 1 will be discussed in further detail later in 
the Results section. The paper now proceeds, though, with a discussion of the theoretical underpinning of 
the analysis presented - an approach by which the Helmholtz viscous vortical flow equation is reduced to 
the classic heat conduction equation. 


II. A New Class of Flow Problem 

A new class of flow problem has been identified that not only encompasses the “moving boundary” 
class of problems (including the classic Lamb-Oseen laminar two-dimensional vortex) but the vortex 
reconnection problem as well. Key to this new class of flow problem is the assumption that the axial and 
radial velocity components, v z and v r , are not merely equal to zero (as in the case of the Lamb-Oseen 
vortex) but the radial velocity is instead proportional to the tangential velocity gradient, e.g. 
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Which dictates, in turn, that 
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For such a class of flow problems (subject to Eqs. 1-2) the convective acceleration terms cancel out the 
vortex stretching terms and thereby achieving the objective of reducing the Helmholtz equation for laminar 
vortical flow (Eq. 3) to the unsteady heat conduction equation (Eq. 4). 
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Note that in the above the term Dcn/Dt is the particle derivative for the vorticity. And, further, the 
vorticity components can be defined by 


co = curFV = V x V 


Subject to the following constraint 


diva) = V • co = 0 
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The result is the employment of length-scale factors to solve for the various three-dimensional 
components of the vorticity and velocity. These length-scale factors are dependent on whether or not the 
inherent functionality (in terms of source functions) is axially symmetrical or non- symmetric. Requiring 
that 
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The following formulation of length-scales holds. Let the axial and radial vorticity analytic functions 
be expressed as a summation of N source terms arranged in spatial distribution. 
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Next, assume that a length-scale can be applied to each individual source term such that 
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Making the above substitution, the following demonstrates that the nonlinear terms in the laminar 
axisymmetric version of the unsteady Helmholtz equation cancel out. First, note that 
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Then, if one assumes that a single common length-scale (for the non-symmetric case) can be applied, where 
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It can be readily seen from the above that 
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The above (Eqs. 7-12), therefore, validates the utility of vorticity/velocity length-scale factors to 
reduce the Helmholtz equation to the heat conduction equation. In particular, it demonstrates that 
superposition of elemental functions representing the vorticity is possible as long as the length scale factor 
is uniformly applied to all elemental (source) terms. Therefore, this also validates the intended use of semi- 
continuous distributions of point and finite- core source terms to model the cut, or broken, vortex filaments 
in the reconnection problem. 

Note, though, that when spanwise symmetry is preserved (as is the case for the vortex filament with 
two breakpoints in close proximity to each other), then there are very special implications as to the length- 
scales employed for, and applied to the analytic functions describing, the vorticity distribution). Assume, 
first, that the equation describing the vorticity distribution can be decomposed into two elemental functions, 
such that 
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Where these two elemental functions can be defined (to the detail required) as follows: 
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Further the vorticity symmetry condition 
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Given axial flow symmetry considerations 
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And assuming that the two length-scales (£ A and i B ) need be applied to relate velocity to vorticity such 
that 


This, in turn, dictates that 
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Canceling out the vorticity term on both sides of the equation yields the length scale constraints for 
symmetrical flow 
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The above (Eqs. 13-20) holds for spanwise symmetrical flow where the plane of symmetry is located 
at z=0. In the above, each (sign of the) length-scale can only be applied to the associated analytic vorticity 
function on the same side of the plane of symmetry in which the length-scale is defined. In the context of 
the vortex filament reconnection problem, the non- symmetric flow length-scale constraint applies to the 
vortex breakdown for a filament with isolated breakpoint or an unbroken vortex filament; the symmetric 
flow length-scale constraint applies to the vortex filament reconnection process where the two breakpoints 
are relatively close proximity. The symmetric flow length-scale constraint also allows for the possibility of 
multiple vortex breakpoints with varying segment lengths, as along as spanwise symmetry is preserved. 
Finally, the above length scale relationship, though derived on the basis of axial velocity and vorticity 
symmetry, is also applicable to definition of the radial velocity. I.e., the following (Eq. 21) must also hold 
(for problems with symmetrical axial velocity and vorticity) 
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Consequently, for this class of flow problems, where length scale factors can be applied, the unsteady 
viscous Helmholtz vorticity equation reduces to the well-known unsteady heat conduction equation. This 
analytical reduction of the vorticity equation to the unsteady heat conduction equation was previously 
demonstrated in Ref. 12. There is a large body of work related to the analytical solution of the heat 
conduction equation, including the use of unsteady finite- volume source modeling, that enabled the 
derivation of (albeit laminar) solutions for the vortex reconnection problem. 


III. Basic Framework for Analytical Solutions 

Given the above introductory discussion, the general solution approach for the vortex reconnection 
problem, assuming axisymmetric flow, is now outlined. This same solution approach was applied to all 
four initial vorticity distributions studied in this paper. Following this outline of the basic analytical 
framework, the four different vortex source terms used in the modeling of the vortex reconnection problem 
is presented. Because of space limitations only the full solution for the instantaneous vortex filament 
reconnection problem will be presented in the Appendix. However, as will be readily seen, the work 
presented in the Appendix is of general applicability to the other three initial vorticity distributions. The 
details of these other analytical solutions can be found in Ref. 12. 

1. The axial vorticity distribution can be described by a semi-infinite continuous distribution of 
unsteady source terms: 
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2. In turn, using length scale factors, the axial velocity can be determined from the axial vorticity; 
this is, in part, dependent on whether the flow is considered nonsymmetrical or symmetrical with respect to 
the filament axis. 


Nonsymmetrical (Eqs. 9b and 11): 
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Symmetrical (Eqs. 17 and 20): 
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3. Given the above derivation of the axial velocity, the tangential velocity can be derived by means 
of the linear ordinary differential equation (derived from Eq. 5). 
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The radial vorticity can now be derived (based, again, on Eq. 5). 
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If the flow is symmetrical, the radial vorticity expression will need to be further separated into the two 
terms, co rA and co rB , such that 
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5. Next, the radial velocity, depending on whether the flow is considered symmetrical or 
nonsymmetrical, can be derived. 

Nonsymmetrical (Eqs. 9a and 11): 


v r = - £co r 


Symmetrical (Eqs. 20 and 21): 
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5. Finally, the tangential vorticity can be derived to complete the solution for the unsteady, three- 
dimensional vorticity/velocity flow components (derived from Eq. 5). 
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As noted above the solution of the vortex reconnection problem is dependent upon the type of source 
term used in the heat conduction equation. In this section four different source terms are briefly 
summarized that correspond to the four different initial vorticity distributions studied in this paper. 


A. Instantaneous Vortex Filaments (Point Source Term) 

Reference 1 6 provides a wealth of information regarding the solution of the heat conduction equation, 
including the use of point and line sources to solve for unsteady problems. Among the source terms 
identified in Ref. 16 is the well-known point source. Equation 22 employs point sources, Eq. 29, to define 
semi-infmite “instantaneous” vortex filament segment expressions for axial vorticity. 
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The complete vortex reconnection solution for instantaneous vortex filament segments is summarized in 
the Appendix. 


B. (Slightly) Parabolic Finite-Core Vortex Filaments 

Reference 16 provided exact expressions for spherical finite-volume source terms; one such expression 
is noted below 



(31) 


The parameter r Q in Eqs. 30 and 31 is a spherical finite-volume source term initial “effective” radii. 
For the above source term, Eq. 30 or Eq. 31, r Q is well defined, i.e. co z \ = a nonzero constant for 

'uOUVCC' 

r <r Q and co z \ = 0 for r>r Q . Correspondingly, r c0 is an initial “finite core” vortex filament core 
radii, such that 8vq /dr\ = 0 for t = 0. In general, for finite-core vortices, with the notable exception of 


= r c0 


uniform-core vortices, the underlying source terms’ initial “effective” radii is not equivalent to the initial 
vortex core size, i.e., in general, r c0 * r Q . There is no a priori relationship between r Q and r c0 , except for 
the uniform-core (constant vorticity in the core) case where r Q = r c0 . As it is anticipated that resulting 
parabolic-core vortex filament vorticity distribution will be nearly uniform, it is assumed - with the 
predictions in this paper made accordingly - that r Q * r c0 . The validity of the resulting parabolic-core 
vortex filament analytic solution (as detailed in Ref. 12) is unaffected by this assumption other than a small 
effect on the spatial-temporal scaling of the predicted flow phenomena, e.g. as will be seen later the core 
size growth trend for the parabolic core solution is slightly under-predicted with respect to nondimens ional 
time. 



C. Dual-Core Vortex Filaments 


Reference 16 also introduced an “approximate” finite volume point source expression (Eq. 32). This 
finite volume source is approximation to the spherical finite volume source term noted in Eq. 31, valid for 
small values of r Q . This “approximate” source term is, as demonstrated in Ref. 12, an exact particular 
solution of the unsteady heat conduction equation. Therefore, though in one sense an approximation, it is 
perfectly valid to employ the Eq. 32 source term to define a “dual-core” vortex filament solution wherein 
the vortex has inner core of negative vorticity. Not unexpectedly, this dual-core solution has a lower limit, 
time-wise, as to its potential physical validity as singularities in the vorticity/circulation expressions are 
manifested as t = 0 . This dual-core model does exhibit predicted vortical flow characteristics that have 
analogous properties to real flow phenomena (e.g. variable vortex core axial distribution and the 
manifestation of vortex breakdown bubbles) and as such is included in this paper. 



(32) 

Note that the parameter r 0 is still explicitly incorporated in the Eq. 32 source term expression because 
of its underlying heritage with respect to the Eq. 31 source term (of which it is a small core radius 
approximation thereof); though, as noted above, this parameter is indeterminate at its nominal defining 
initial condition because of the manifestation of a singularity in the source term at t = 0. However, given 
the earlier proposed initial condition approximation, i.e. r Q « r c0 , for the parabolic core vortex filament 
source term, an equivalent parameter value “assignment” logically can be made for the dual-core source 
term, i.e. r Q => r c0 ; by this means it is asserted that physically meaningful predictions can be made at all 
but the smallest values of nondimensional time. 


D. Uniform Core Vortex Filaments 

The following source term yields a solution wherein the vorticity radial distribution can be defined as 
initially uniform (versus the slightly parabolic initial distribution of the solution derived using spherical 
finite-volume sources). From Ref. 16, for the placement of the point sources at arbitrary polar-cylindrical 
coordinates (versus at r=0 for the “instantaneous” vortex filament segments), the following holds 
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The above source term can be employed to define semi- infinite continuous line sources, subject to 
integration with respect to 6* (from 0 to 2jt) and z* (piecewise from -s to -oo and s to oo), that then can 
be integrated with respect to r * (from 0 to r c0 ) to yield a uniform-core vortex filament solution. Equation 
33 integration with respect to 6* closely follows the work detailed in Ref. 16 and yields a modified Bessel 
function of the first kind of zero th order, 7 0 (x) . Correspondingly, integration with respect to z* yields 
results similar to the instantaneous vortex filament solution, which is detailed in the Appendix. The 
uniform-core solution is unaffected by the earlier assumed relationship, or rather equivalence, between r Q 
and r c Q . Complete solution details for the uniform-core vortex filament solution are provided in Ref. 12. 



The resulting uniform-core solution has been found to have excellent agreement with an alternate related 
analytic treatment detailed in Ref. 17. 


IV. Discussion of Results 


A. Time-Dependent Behavior of Unbroken Vortices 

In addition to the well-known Lamb-Oseen vortex model, there are a number of alternate unsteady two- 
dimensional vortex models noted in the literature - e.g. Refs. 12, 17, 18. Figure 6a-b compares the 
tangential velocity profiles, as a function of time, of the uniform finite-core solution (by way of illustration) 
to alternate unsteady vortex models cited in the literature. For example, the uniform finite-core solution is 
compared to both the Lamb-Oseen and Kirde (for two different parameter values of n = 1/2 and n = 3/4) 
vortex models. The Kirde solution (Ref. 18) is based, in part, on confluent hypergeometric functions. 
Figure 6a-b also compares the uniform finite core solution to the Aboelkassem, et al, unsteady, two- 
dimensional Rankine-like vortex model (Ref. 17). (It should also be noted that the Ref. 17 methodology 
was later extended to other types of “monopolar” vortices, see Refs. 19-21.) 

The Aboelkassem model and this paper’s parabolic and uniform finite-core vortex solutions all 
transition from a Rankine-like initial tangential velocity profile to a Lamb-Oseen profile over time. The 
Aboelkassem vortex model is based on a Fourier-Bessel series formulation. The parameter values 
n = 1000 (number of Bessel series terms) and /5 = 1/100 (a multiplier constant) are used in the 
Aboelkassem model, as recommended. As the uniform finite core solution and the Aboelkassem vortex 
model both implicitly assume a Rankine vortex as an initial velocity profile, it is anticipated that there 
should be very good correlation between the two models, irrespective of their mathematical 
treatment/formulation. This is indeed the case. For all but the very earliest of nondimensional times, the 
uniform finite-core solution and the Aboelkassem models show excellent agreement. Both solutions have 
their strengths and weaknesses, though. The advantages of the Aboelkassem model is that it does not 
exhibit any numerical accuracy or stability issues throughout the complete range of solution applicability 

( 0 < t < oo). The main disadvantage of the Aboelkassem model is that a non-rigorous “calibration” of the 
Aboelkassem vortex model with the Lamb-Oseen velocity profile was performed (at a single point-wise 
nondimensional time, which, in Ref. 17, was somewhat arbitrarily selected) in order to define the 
recommended/suggested value for the parameter /3 . The uniform finite-core solution doesn’t have to resort 
to such a mathematical artifice in its derivation. However, the uniform finite-core solution does exhibit 
numerical stability issues for very small values of time in its current mathematical formulation and 
numerical solution implementation. This can be partially compensated for by increasing the order of series 
terms included in the solution; the more proactive solution to the numerical stability issue would be to 

derive a specialized asymptotic formula for the uniform finite-core model for t* -> 0. 
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Fig. 6. Comparison of Unsteady Vortex Models: (a) t = 0.01 and (b) t =0.1 
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The core radius growth (relative to the initial core size) with time for the various vortex models can be 
seen more clearly in Fig. 7 a. The three vortex models (the Kirde (n = 1/2 and n = 3/4), Aboelkassem, and 

uniform finite-core models) exhibit core radius growth with time as approximately r c 'yft . The 
Aboelkassem and uniform finite core models’ core radius growth trends are reasonably close in agreement 

with the modified Lamb-Oseen core radius growth trend of r c = ^r c 2 0 + 4 avt where a = 1.25643 is the 
Lamb constant. 


The uniform and parabolic finite-core tangential velocity profiles are nearly identical as discussed 
earlier (except at the earliest nondimensional times). The influence of the dual-core, parabolic, and 
uniform core models, with respect to core growth with time, is shown in Fig. 7b. As can also be seen in 
Fig. 7b both the parabolic and dual-core vortex core size growth trends are under-predicted (the parabolic 
trend more so then the dual-core). This is likely a consequence of the “near-uniform-vorticity” assumption 
used in defining an “equivalence” between r Q and r c0 . This equivalence assumption will have to be 
reassessed in future work. Nonetheless, this issue does not fundamentally affect the overall validity of the 
parabolic and dual-core predictions (and underlying analytical solutions), it is instead primarily a 
spatiotemporal scaling issue and not a question of the basic flow phenomenology. 
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Fig. 7. Core Radius Growth Due to Viscous Diffusion: (a) Comparison to Models in Literature and 
(b) Comparison of Models Summarized in this Paper 


B. Reconnection and Breakdown 


Figure 8a-b shows the circulation distribution, along the filament axis, as a function of nondimensional 
time. The diffusion-like nature of the reconnection process is readily apparent in this figure as the vortex 
filament segments reconnect in the intermediate region between the vortex breakpoints. Secondary features 
in the circulation, though, can be seen resulting from the inner core of negative vorticity inherent in the 
dual-core initial vorticity distribution. 



Fig. 8. Vortex Circulation Axial Distribution as a Function of Nondimensional Time: (a) 
“Instantaneous” and (b) Dual-Core Vortex Models 


The secondary circulation distribution features observed in the Fig. 8b culminate in singularities (though 

not shown) at the vortex filament breakpoints as t* — > 0. Therefore, a lower bound as to nondimensional 
time is required to prevent prediction of physically unrealistic flow. This lower bound was derived, in Ref. 

12, to be approximately t* > |r* 0 J jlO . 

Enstrophy is the mean square of all three vorticity components of the flow. Enstrophy is estimated by 
the expression f] -> + co^ + cdq J^3. In this particular set of flow problems, axial vorticity is the 

principal component for the enstrophy contours in Fig. 9a-c. The full complexity of the predicted vortex 
reconnection process can be seen in this figure. Presenting enstrophy contours for a sequential set of 

nondimensional times (t* = 0.03 (parabolic- core only), 0.12, 0.2, 0.3, and 1.0) for one vortex-core-size-to- 
cut ratio ( r* 0 = 1.5 ) provides a powerful means of visualizing the vortex reconnection process. At the 
earliest nondimensional time values, the tangential and radial vorticity components contribute more 
significantly then at later time-wise. It is also at these early moments of the vortex reconnection process 
that clear differences can be distinguished between the enstrophy contour predictions of the three finite- 
core initial vorticity distributions. At later times, the enstrophy contours become nearly indistinguishable 
from each other. The axial flow wave front propagation boundaries can also be clearly seen in Fig. 9a-c as 
abrupt discontinuities in the enstrophy contours. 



Fig. 9. Enstrophy ( r c0 =1.5, Re=50, equivalent scaling, and at, sequentially, time values of t = 0.03 

(Parabolic-Core Only), t = 0.12 , t = 0.2 , t = 0.3 , and t = 1.0)): (a) dual-core solution, (b) 

parabolic-core, and (c) uniform-core 


Figure 10 is a typical predicted velocity vector map (axial and radial velocities) superimposed over a 
color contour plot of the tangential (or swirl) velocity distribution in the intermediate region between the 




(two) “cut” vortex filament breakpoints. This velocity flow field, in the intermediate region between vortex 
filament breakpoints, is typical for all initial vorticity distributions studied. 



Fig. 10. Vector Field (vector components are v r and v z ; v 0 indicated by color/shaded contours) for the 
r/z Plane (t* = 0.1) 


A universal flow feature of the vortex reconnection process appears to be the formation of the “funnel- 
like” stream function contours and spiral streakline, or particle trace, flow patterns. Figure 11 presents a 
composite image of both the stream function contour and particle trace patterns for a typical result. 
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Fig. 11. Composite Image of Stream Function Iso-Surfaces (t* = l and Re=50) and Spiral Flow 
Pattern (Particle Trace between 1 < t <100) 



Both the parabolic and the dual-core initial vorticity distributions yielded predictions of variable core 
size distributions with respect to the axial coordinate (Figs. 12a-b and 13a-b). The dual-core predictions 
show greater variation in the core radius than does the parabolic initial vorticity distributions. It is 
interesting to observe that variable core radius distribution is not, as might be initially expected, solely a 
consequence of reduced or negative vorticity near the center of the vortex core, as the parabolic initial 
vorticity distribution also exhibits this flow feature. The under-prediction of the dual-core and parabolic 

core size in the far- field ( z* -> oo) for these relatively small nondimensional times can also be seen. 
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Fig. 12 . Variable Vortex Core Size: (a) t* = 0.12 and (b) t* = 0.2 
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Fig. 13. Three Dimensional Contours of Variable Core Size: (a) ( t = 0.12), and (b) ( t = 0.2) 


As noted in Table 1, for the dual-core initial vorticity distribution, the vortex filament segments 
manifests flow features similar to that observed for vortex breakdown phenomena including the sudden 
expansion of the vortex core and the formation of recirculation bubbles. Figure 14a-d presents typical 
prediction results for the vortex core radius and breakdown bubble envelopes. Figure 15a-b goes further by 
illustrating the predicted vortex breakdown bubble geometry trend (as represented by both the bubble 
fineness ratio and the axial location, as a fraction of bubble length, of the maximum bubble thickness) as a 
function of nondimensional time and initial vortex core size. The predicted vortex breakdown bubble goes 
from being squat, acorn-like, to slender nearly spindle-like. Similar bubble topologies have been 
experimentally observed. Unique, perhaps, to the predicted vortex reconnection process, though, is that 
these vortex breakdown bubbles only occur in the early stages of the reconnection process and then 
dissipate/disappear after a critical nondimensional time value is reached; thereafter, the reconnection 
process is insensitive to the initial vorticity distribution and proceeds, in general, in a fairly similar manner. 




Fig. 15. Bubble Geometry Trend as a function of nondimensional time: (a) fineness ratio and (b) 
axial location (as a fraction of bubble length) of maximum thickness 


Given the mathematic character of the dual- core vortex filament flow equations, it is relatively 
straightforward to identify discrete axisymmetric vortex breakdown bubble “finite states,” and the 
nondimensional time limits that conform to them, that are exhibited in the predictions. Considerable 
discussion, including qualitative comparisons with experimental measurements and observations in the 
literature, is given in Ref. 12 as regards the prediction of this analogous flow behavior to the classic vortex 
breakdown phenomenon. 


V. Concluding Remarks 

Analytical solutions for the initial breakdown and subsequent reconnection of a laminar vortex filament 
- for the ideal case of having been “cut” or “broken” subsequent to a collision with a solid body — has been 
developed. The analytical solutions do not treat the actual collision process itself, though. These analytical 
solutions belong to a general class of flow problems, identified through this work, that encompasses as a 
subset the extensively studied “moving boundary” laminar flow problem. The viscous Helmholtz equation 
for vorticity is reduced to the unsteady heat conduction equation though a novel treatment of the problem. 
Consequently, instantaneous and/or finite-core line sources for the vorticity can be used to model the 
unsteady flow vortex filament breakdown and reconnection processes. A diverse and fascinating catalog 
of flow features is predicted by the resulting analytical solution - including several examples of analog 
vortex breakdown flow structures (bubble and conical types) similar to that previously reported in the 
literature. The breadth of flow phenomenology predicted by the four initial vorticity distributions examined 
in this preliminary work begs continued investigation in this area, including a comprehensive examination 
of the influence of initial vorticity distributions on vortex evolution. 

It is anticipated that the work presented in this paper will find future application to refined rotary- and 
fixed-wing computational wake models. Additionally, studies into active control/dissipation of vortical 
structures - such as trying to enhance the decay of trailed wing-tip vortices of large commercial aircraft to 
minimize the occurrence of wake upset for such aircraft - will also hopefully be furthered by this work. 
Finally, the general study of the vortex breakdown phenomena will also benefit from this work. 
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Appendix - Instantaneous Vortex Filament Solution 

Reference 16 reveals that a particular solution of Eq. 4, for an instantaneous point source, is given by 
(Eq. 29 presented again for clarity) 


co 


z 


point source 



+ z-z* 



Where c is a constant, z* is the origin of the instantaneous point source with respect to the z-axis, and r 
and z are the coordinates of the point for which the vorticity is predicted. 

Correspondingly, Ref. 16 derives the solution to the unsteady heat conduction equation for an 
instantaneous line source of infinite length (represented as a continuous distribution of instantaneous point 
sources) by Eq. 34. 
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infinite line source 
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(34) 

With the constant, c, set equal to the circulation strength, y, (Ref. 15), then Eq. 34 becomes the 
vorticity distribution for a Lamb-Oseen vortex. 

To solve the unsteady laminar flow problem for a rectilinear vortex filament that has been “cut” or 
“broken,” however, it is necessary to represent the vortex filament segments with semi-infinite distributions 
of instantaneous point sources (Eq. 35a-c). Note that the constant, c, has been assigned the value c=y. 



(35a-c) 

The integral terms in Eq. 35 can be found in standard handbooks for integration formulas (see Ref. 22). 
The solution of vorticity component with respect to the z-axis is given by 
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Where erfc^x) is the complementary error function, i.e., 



(37) 


Noting that complementary error function is related to the error function by the relationship 
erfc{x ) = 1 - erf(x) . And, further, erf(-x ) = -erf{x) , erf(o o) = 1 ? an d erf[0) = 0 . With these properties of 
the error function and complementary error function in mind, it can be readily seen that Eq. 36 satisfies 
three important boundary constraints and/or conditions. 

a) ^-^— e ~ r2 ! Avt For (a) s -» 0 

A jtvt 


Or when 



Or when 


(c) j 0 And |z| ^ ^ 


(38a-c) 

The vorticity distribution approaches that of the Lamb-Oseen vortex as the above conditions (Eq. 38a- 
c) are met or approached. 

Having derived the vorticity distribution, the axial velocity distribution follows straightforwardly from 
Eq. 23b, given Eqs. 35-36. 


Where 
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The appropriate length scale factor in Eq. 39a-c is the separation distance, specifically, £oc-s. Besides 
being a natural length scale factor, using s as the length scale automatically satisfies the “moving 
boundary” condition of the vortex filament breakdown/reconnection problem; i.e., when s = 0 then v z = 0 . 
This choice for the length scale also satisfies that constraint that v z > 0 for z > 0 and v z < 0 for z < 0. 
However, the vortex filament solution also has to take into account a vortex axial flow “wave front” 
propagation boundary moving at a finite velocity, V P , along the axis of the filament segments. Therefore, if 
\z\^ s + Vpt then v z ^ 0 ; if \z\ > s + V P t then v z = 0 . It is especially important to note, though, that for the 

instantaneous vortex filament solution, this vortex axial flow wave front propagation occurs effectively 
instantaneously, i.e. V P -> oo ? such that the axial velocity equation does not have to account for a wave 

front boundary cut-off.) Therefore, in general, the length scale now becomes £ = 
where is the unit (Heaviside) step function. An inherent assumption is that the “wave front” 

boundary axial location varies with time but not the radial coordinate. In general, the following definition 
for the wave front propagation velocity is employed 



(40) 

Equation 40 assumes that the axial flow established by the vortex filament cut/break is primarily contained 
within the vortex core (which is only an approximation). Additionally, note that the integration of the axial 
velocity profile occurs behind and not ahead of the wave front. Specifically, though, noting the definition 



of the wave front propagation velocity, Eq. 40, that for the instantaneous vortex filament case V P — > oo 
and, therefore, £ = -s- 

Application of the length scale factor to the axial vorticity expressions gives the following expression 
for the axial velocity 
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Given Eqs. 24 and 36, the tangential velocity distribution can be derived through solution of a first- 
order ordinary differential equation. 
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Which has the general solution (Ref. 22) 
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Performing the required integration yields the expression 
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The constant C (with respect to r) can be found by imposing the constraint that v 6 = 0 for r=0, for all z 
and t. Therefore, 
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And so the final solution for the vortex filament tangential velocity distribution is 
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The analogous attributes of Eq. 46 with respect to the Lamb-Oseen vortex are obvious (Ref. 15). As 
demonstrated in a similar manner with respect to the vorticity distribution, if ^ — > 0 then, correspondingly, 
the tangential velocity distribution reduces to the Lamb-Oseen profile. Further, note that the vortex core 
radius is equivalent to that of a Lamb-Oseen vortex, i.e. r c = 2.24 ^fvt , and is invariant with respect to the 
filament break, in terms of s, or the axial coordinate, z, along the filament segments. 

Next, given Eqs. 25, 27b, and 46, the radial vorticity, co r , and velocity, v r , can be derived by taking the 
derivative of v 0 with respect to z. Noting first, though, that the following holds true 
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Therefore, applying the above formula for the derivative of the complementary error function the 
following results 
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Applying the length scale factor, i = -s, to Eq. 48, given Eq. 27b (presented below, again for clarity), 
such that 


v r ~ ^ A^r A ^ B^rB 
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This gives the following expression for the radial velocity distribution 
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A key outcome of Eq. 49 is the prediction that the radial velocity is negative (inward to the vortex 
core) for all values of r, z, and t. The radial inward flow contributes to the manifestation of “funnel- like” 
stream function contours and spiral (inward from the vortex breakpoints) streakline, or particle trace, flow 
patterns. This general vortex reconnection flow behavior holds true for all of the initial vorticity 
distributions studied in this paper. 



